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Abstract. Predictions for the total cross section for e+e^ — > tt near threshold are reviewed. The 
renormalization group improved results at NNLL order have improved convergence and reduced 
scale dependence relative to fixed order results at NNLO. Prospects for measurements of the top- 
quark mass, width, and Yukawa coupling are discussed. 



INTRODUCTION 



O . The production of top quark pairs is among the important projects of a future linear 

collider. The large top quark width ~ 1.5 GeV makes the threshold cross section look 
quite different from that of charm or bottom pairs. The large value of the width prohibits 
O . the production of toponium states, and at the same time serves as an infrared cutoff 

from sensitivity to non-perturbative effects. Thus, perturbative methods can be used to 
. ^ describe the top-antitop dynamics to a very high degree of precision. This makes the 

p ■ threshold region an ideal place for extracting fundamental top quark parameters such as 

the top mass, width, and Yukawa coupling (for a light higgs). 

In the threshold region ^/s ^ Inif ± 10 GeV, the t and t move with non-relativistic 
^ ■ velocities. Defining the energy m,v^ = ^/s — lnit we see that this region of s corresponds 

to velocities |v| ^ a^. In this region an exact treatment of QCD Coulomb singularities 
\ [cLs/^Y is required, ruling out a pure a? expansion. However, a combined expansion in 

5^ ■ powers of v and can be performed. Schematically, for the e^e^ tt cross section, 

(3tt{s), we want an expansion of the form 

k 

(1) 



= ^ = v£(^) X [l + {a.,v} + {a?,a.v,v2} + ... 



■'^J^/' k 



LO NLO NNLO 



The power counting for these corrections can be implemented in a simple and systematic 
way using the effective theory framework of Non-Relativistic QCD (NRQCD). 

The leading order (LO) prediction for R is shown in Fig. 1. The large top width gives 
the threshold cross section a smooth line-shape, with a single bump from the remnant 
of the IS toponium state. The characteristics of the cross section are sensitive to the 
top parameters. In particular, the top mass determines the location of the rise/peak. 
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FIGURE 1. Leading order predictions for R. Results are shown for three renormaHzation scales of order 
the momentum transfer in the Coulomb potential, /j = 15 GeV (upper), 30 GeV (middle), and 60 GeV 
(lower). 




the top width determines the slope of the rise and shape of the peak, and the overall 
normalization provides information on as{mt) and the top- Yukawa coupling for a light 
higgs. However, at LO the prediction for R suffers from considerable scale uncertainty 
as shown by the three curves in Fig. 1. This illustrates the importance of including 
higher order terms. Several groups have computed the next-to-next-to-leading order 
(NNLO) QCD corrections to the total cross section (see Ref. [1] for a summary and 
comparison of these results). Using "threshold" top-quark mass parameters, an accurate 
prediction for the location of the peak in the cross section was obtained. It was concluded 
that infrared safe top masses could be determined with a precision of < 200 MeV. 
However, surprisingly the NNLO corrections to the cross section normalization were 
big, with large scale uncertainty. The residual scale dependence was estimated to make 
the normalization of the cross section uncertain to ~ 20% [1]. These large corrections 
seemed to jeopardize measurements of the top width, strong coupling, and Yukawa 
coupling. 

To understand the source of this large scale dependence it is useful to recall that the 
dynamics of the top-antitop system are governed by vastly different energy scales. These 
quarks have mass ~ 175 GeV, typical momenta p ~ /%v ~ 25 GeV, and energies 
: m^v^ ~ 4 GeV. The scattering amplitudes therefore involve logarithms 



'"(S)' '"(S- 



which are not all small for a single choice of /j. In the expansion for R these large 
logarithms appear as logs of the velocity. This suggests that a better expansion might 
involve a summation of these logarithms 



^ = )^£(a,lnv)'x [l + {a,,v} + {ala,w,y^} + . . . 



k 



(3) 
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The additional summation of logarithms can be performed using renormalization group 
equations (RGB). This is similar to how logarithms of m-w /mjj are summed for the elec- 



troweak Hamiltonian. In our case the summation is complicated by the presence of two 
low energy scales that are coupled by the equations of motion, E — p^/ (2mf). This com- 
plication can be dealt with by using a renormalization group with a subtraction velocity 
V [2, 3, 4, 5, 6]. In this framework there are two renormalization group parameters in 
the effective Lagrangian: i^s for soft gluons and idy for ultrasoft gluons, where jds = mtV 
and fuj = nitV^. Running from v = 1 to V ~ v sums logarithms of v and minimizes both 
ln(/j|/p^) and ln(/j^/£'^) terms in the amplitudes. (The ln{ij^/mf) terms are minimized 
by matching QCD onto the effective theory near /j = nit.) In Refs. [7, 8] it was shown 
that the expansion in Eq. (3) is better behaved than that in Eq. (1). Furthermore, the 
normalization of the cross section at NNLL has significantly smaller scale uncertainty 
3%) than the NNLO result. This review focuses on these results; for more details and 
further references see Ref. [8]. 

Below, the steps in the renormalization group improved calculation are briefly de- 
scribed. This is followed by predictions for the cross section as well as a description of 
the prospects for measurements of the top parameters in light of the reduced theoretical 
uncertainties at NNLL order. 



RENORMALIZATION GROUP IMPROVED CALCULATION 

The computation of the renormalization group improved cross section can be divided 
into three parts: 

1. Matching of QCD onto an effective theory for non-relativistic top quarks. Deter- 
mine the Wilson coefficients C(v) at v = 1 (/j = nit) to the desired order in as{nit). 

2. Scaling C(v) from v = 1 to v = vq ^ CfUs by calculating anomalous dimensions 
and using the renormalization group. This scaling sums the terms of the form 
[a,ln(v)]^'inEq. (3). 

3. Computing the cross section using the Lagrangian and currents renormalized at the 

low scale v = vq . 

Each of these parts will be discussed in turn. 

1. Matching onto the Effective Theory 

For non-relativistic scattering the relevant momentum regions can be classified by 
their typical energy and momenta (A:°,k). They include hard modes with momenta 
~ (m,m), potential modes with momenta ~ (mv^,mv), soft modes with momenta ~ 
(mv,mv), and ultrasoft modes with momenta ~ (mv^,mv^). Fluctuations involving hard 
or offshell momenta are integrated out, while effective theory fields are introduced 
for modes with nearly on-shell momenta. The degrees of freedom therefore include 
potential top and anti-top quarks (\|/p and Xp), soft gluons and light quarks (Aq and (p^), 
and ultrasoft gluons and light quarks (A^ and Soft energies and momenta appear 
as labels on the fields while ultrasoft momenta are represented by explicit coordinate 



dependence [2]. This enables us to distinguish the size of momenta, for instance a 
derivative 3^\|/p(jc) ~ m\^\\fp{x). 

In this framework the action for non-relativistic top quarks has terms 



L 



2m 



¥pW + (¥^X) 



- L^(p,p') [vI/WvpWxIp'Wx-pW 



p,p 



-27ca,(mv) £ ¥p' [^J,^?]f/ap¥p + (¥ ^ X) , 



(4) 



P,P',9,9' 



where color and spin indices are suppressed. The covariant derivatives in the first line 
involve only ultrasoft gluons. The function F(p, p') in the second line contributes to the 
potential between quarks and anti-quarks. For our purposes 
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(P^ + P'') 



^(P,P') = 7jM^) + + ^'^2{V) + ^1^s{v) + 



mlkl nfi nfi Irn^Vi?' 



'^^r(v), (5) 



where k = p' — p, S is the total spin operator, and all coefficients are in the color singlet 
channel. The matching for these Wilson coefficients is needed at two loops for i^^^ 
one loop for V]^, and tree level for 1^2,5,r [5, 6]. Finally, the third line in Eq. (4) is an 
example of the type of interaction that occurs between potential quarks and soft gluons, 
with f/ap(P:P')?:?0 ^ matching function of the label momenta. At NNLO (and NNLL 
order) time ordered products of two of these soft interactions also contribute terms to 
the potential giving 



Vsoft(P,P') = 



(6) 



where P, are coefficients of the QCD |3-function and the constants a; can be found in 
Ref. [9]. The complete potential is then V(p,p') = F(p,p') +ysoft(p,pO- 

We also need to take into account that the top quarks decay. We will assume that 
the decay products are hard and can be integrated out. At lowest order this induces the 
operators 



(7) 



In the Standard Model the dominant decay channel h t ^ bW^ and gives a width of 
Tt = 1.43 GeV which we will use as our central value. Counting Ff ~ m^v^, Eq. (7) 
gives a consistent next-to-leading order treatment of electroweak effects for the total 
cross section [10, 1 1]. Thus, we will not include electroweak decay related effects to the 
same order as the QCD corrections. From partial knowledge of these corrections [12], 
the missing terms are expected to be at the few percent level. 



Besides the effective Lagrangian we also need external currents to produce the top- 
antitop pair. Since we wish to describe e+e — > {7*,Z*} — > ?f these currents are induced 
by both electromagnetic and weak interactions. The relevant vector current is Jp = 
ci(v) Op,i +C2(v) Op 2, where 

Op,i = \Vp^c{ic2)x%, Op,2 = ^VpVi<y2)X*-p, 
and the relevant axial- vector current is Jp = C3(v) Op^3, where 

Op,3 = ^¥p^[o,op](i02)X-p. (8) 

The matching for these Wilson coefficients is needed at two loops for ci and tree level 

for C2.3. The two-loop matching for ci is scheme dependent [13, 14, 15], and in the MS 
scheme with our definition of the operators, can be found in Ref. [8]. 



2. Renormalization group scaling 

To sum the (a^lnv)-^ terms in R we must determine the anomalous dimensions for the 
Wilson coefficients 1^c,k,2,s,r in Eq. (5) and the current coefficients 01,2,3- The anomalous 
dimensions for l/c^) '^h and Vi^s.r are required at three, two, and one loop respectively. 
These have been computed in Refs. [3, 5, 6, 16, 17], and due to mixing depend on the 
one-loop running of the HQET terms up to 1/m^ [18]. These anomalous dimensions can 
contain terms like (with bi coefficients that depend on color factors) 



- bia^{mtv) + b2asimtv)as{mtv^)+b3a^{mtv)ln — - — - 

+Z74 4(v)a2(m;V) + .... (9) 

For l^c,k,2,r both soft and usoft loops contribute, so their anomalous dimensions depend 
on both as (mfV) and as (mfV^). The anomalous dimension for 'P'v comes only from soft 
loops but also involves the mixing of the \j/o • B\|/ operator whose Wilson coefficient 
is (v). Of the current coefficients, C3 has no anomalous dimension, while C2(v) has 
contributions only from ultrasoft loops [7]. More interesting is the anomalous dimension 
for ci (v) which starts at two-loop order from purely potential loops. At this order [2] 

^'"l^'Ml = -^(^ + ^^(v) + ^^(v)+2n(v))+^if*. (10) 

so the solution for ci depends on the solutions for the potential coefficients. At three 
loops there are new contributions to the anomalous dimension for ci(v) coming from 
mixed potential-ultrasoft and potential- soft loops. These contribute at NNLL order but 
are currently unknown. In Ref. [8] these unknown terms were estimated to affect the 
cross section at the 2% level (this rough estimate was based on the size of known terms 
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FIGURE 2. Running of the color singlet (s) and octet (o) Wilson coefficients for the Jt^/(m|k|) potential 
from Ref. [5]. The solid vertical Une marks the Coulombic region where V ~ vq. 



and using dimensional analysis and parameter dependence to estimate the size of the 
contributions which are unknown). Of all the Wilson coefficients the one which is the 
most responsible for the difference in the NNLO and NNLL cross sections is '^'^(v) 
which changes by an order of magnitude between v = 1 tov = 0.15. The two-loop 
running of this coefficient is shown in Fig. 2. 



3. Cross section computation 

The total cross section for e^e~ tt is given by 



Ana? 

35 



F"{s)R"{s)+F''{s)R''{s) 



(11) 



where and F'^ are trivial functions depending on the charge and weak isospin of the 
fermions, and sinGvj^. and R'^ are determined by 



An 

R\s) = — Im 
s 



R^is) 



An, 



-Im 



Ci(v)jli(v,m,v) + 2ci(v)c2(v) J^2(v,m,v) 



c^(v)jl3(v,m,v) 



(12) 



where Sli are time-ordered products of effective theory currents [q— {y/s — 2mt,0)] 







^2 



(v,m,v) = ^£ [d'xe'^-^(0 



J?3(v,m,v) = iY, d 



Op,iwo;,2(o)+Op,2w o;,i(o) 

7- Op,3(x) 3(0) 







(13) 



These time-ordered products can be evaluated in terms of non-relativistic Greens func- 
tions to give 

J?i(v,m,v) = 18 G^{Y,m,v) + {^2{v)+2^,{v))dG^{Y,m,v) 

+ 't^,(v)5G''(v,m,v) + 't^fc(v)5G^(v,m,v) + 5G'''"(v,m,v) 

Jl2(v,m,v) = v^J^i(v,m,v) , J^3(v,m,v) = 12G^(v,m,v)/m^, (14) 

Here are evaluated numerically with the 1 /k^ term in F(p, p') and Vsoft(p, p') [19]. In 
Ref. [8] we analytically evaluated SG^'*"'^'^'" with a single insertion of the corresponding 
potentials or p"^ kinetic energy correction. The P-wave Greens function G^ was also eval- 
uated in closed form. The analytic calculations enabled all ultraviolet subdivergences 
to be subtracted in MS which is necessary to be consistent with the scheme depen- 
dence of the Wilson coefficients. In Eq. (14) the velocity v = — Init + zTj] '^1'^ jnij'^ 
and m — nit is the pole mass. The Greens functions depend on the subtraction veloc- 
ity through ln(v^/v^) and these logarithms are not large when the Greens functions are 
evaluated at the low scale v ~ vq. At this scale all large logarithms have been resummed 
in the potential and current Wilson coefficients. Typically, vq ~ 0.15 — 0.2, but to nu- 
merically test the remaining scale dependence we use the larger range vq = 0. 1 — 0.4. 



RESULTS 

Soon after the NNLO results were derived it was realized that the inherent uncertainty 
in the top-quark pole mass due to infrared renormalons causes problems for predictions 
for the peak in the cross section. Therefore, for precision predictions the pole mass is 
not a suitable mass parameter. The MS top-mass provides an infrared safe alternative, 
however it complicates the non-relativistic power counting. Essentially, it shifts E by 
an amount ~ m^a^ which is much larger than the original size of the energy £ ~ mya^. 
Both of these problems can be addressed by switching to a "threshold mass", defined as 
an infrared safe mass parameter which differs from mf°^^ by ~ mtO^. Several possible 
threshold masses were suggested, including the PS mass [20], kinetic mass [21], and IS 
mass [22, 12]. In Ref. [1] it was concluded that threshold masses could be determined 
with a precision of < 200 MeV from the total cross section. Converting the result to an 
MS top-mass would then lead to a similar precision for this parameter. In this section 
predictions will be given using the IS mass parameter. For a detailed description of how 
the NNLL pole mass expressions are converted to the IS mass in a manner consistent 
with the power counting see Ref. [8]. 

We begin by comparing results in the fixed order and renormalization group improved 
expansions. We concentrate on since gives only a small contribution to the cross 
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FIGURE 3. Comparison of QjR'' with fixed m} ^ mass for the fixed order and resummed expansions. 
The dotted, dashed, and solid curves in a) are LO, NLO, and NNLO, and in b) are LL, NLL, and NNLL 
order For each order four curves are plotted for v = 0. 1, 0. 125, 0.2, and 0.4. 
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FIGURE 4. (a) Position of the peak in the cross-section versus v at NNLO (dashed) and NNLL (solid). 
The vertical line at v = 0. 15 is a physically motivated endpoint for the running, (b) Variation of the NNLL 
cross section for a ±10% change in the value of the top quark width. 



section, and is essentially identical in the two approaches. The results are shown in 
Fig. 3 and use M}^ = 175 GeV, a,(mz) = 0.118 and = 1.43 GeV. At each order 
in the expansions four curves are shown which correspond to v = 0.1, 0.125, 0.2, and 
0.4. It is clearly visible that the NNLL results in Fig. 3(b) have much smaller scale 
dependence than the NNLO results in Fig. 3(a). It should be noted that our NNLO 
results shown in Fig. 3(a) agree quantitatively with those presented in Ref. [1]. The 
uncertainty in these results stems to a large extent from the uncertainty in the choice of 
the renormalization scales in the NNLO contributions. Essentially what the anomalous 
dimensions and renormalization group do is remove this uncertainty. Also, more than 
half of the improved convergence of the NNLL result is due to the reduced size of 
'P'ki^ = vo) compared to l^k{'^). 

From the remaining scale uncertainty and the size of some higher order QCD cor- 
rections, the uncertainty in the NNLL cross section was conservatively estimated to be 
±3% [8]. This level of precision should enable extractions of various top parameters 
from the cross section with fairly good precision. In Fig. 4 we show the scale depen- 
dence of the peak position for the NNLO (dashed) and NNLL (solid) predictions. The 
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FIGURE 5. Variation of the NNLL cross section for (a) the inclusion of a Standard Model (SM) Higgs 
boson and (b) the value of the strong coupling. The relative changes are shown by red dashed lines. For 
(a) the lower black solid line is the decoupling limit for the Higgs boson, and the upper blue solid line is 
for a SM Higgs with mass mn = 1 15 GeV. 



NNLL prediction is slightly less scale dependent than the NNLO prediction until we get 
to small V. For v < 0.15 the larger scale dependence at NNLL is explained by the fact 
that these predictions depend on the coupling as{mtV^), while the NNLO predictions do 
not. Also shown in Fig. 4 are the NNLL predictions for the total cross section varying 
the width by ±10%. The size of the variations indicate that a measurement with better 
than 10% precision is definitely feasible. In Fig. 5 the dependence of the cross section on 
the top Yukawa yt (for a Higgs mass mn = 115 GeV) and on as(mz) are shown. It looks 
quite promising that a ±20% variation in yj gives a larger change in the cross section 
than our estimate for the remaining theoretical uncertainty. It should be kept in mind that 
since both yt and as{mz) mainly effect the normalization, at some level these parameters 
cannot be fixed independently using only the total cross section. 



CONCLUSION 

In this talk I have discussed predictions for the threshold e^e^ tt cross section at 
NNLO and NNLL order as defined by the expansions in Eq. (1) and Eq. (3). The NNLL 
predictions made in Refs. [7, 8] sum large logarithms of the velocity using results for the 
renormalization group improved Wilson coefficients from Refs. [2, 3, 5, 6]. One missing 
ingredient is the three loop anomalous dimension for c\, for which only partial results 
are known. However, rough estimates indicate that this missing anomalous dimension is 
unlikely to affect the cross section at more than the 2% level [8]. The stability of pre- 
dictions for the peak in the cross section are very similar at NNLO and NNLL, so that 
measurements with 6m/ < 200 MeV for short distance masses are feasible with either 
expansion. The size of the NNLL normalization corrections and variation of the NNLL 
cross section for various choices of the renormalization parameter are an order of mag- 
nitude smaller than the results of earlier NNLO calculations. A conservative estimate of 
the remaining theoretical uncertainty in the total cross section is ±3% [8]. With such 
small uncertainty, measurements of top parameters with uncertainties h<Xs{fnz) ~ 0.002, 
bTf/Tt ~ 5%, and hyt/yt ~ 20% appear feasible. However, realistic simulation stud- 



ies should be done to see how these numbers hold up once effects such as initial state 
radiation, beamstrahlung, and the beam energy spread are taken into account. 
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